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Abstract.  We  conducted  planar  shock  wave  experiments  to  measure  the 
dynamic  loading  of  dry  and  water-saturated  soils  provided  by  the  Army  Re¬ 
search  Laboratory  (VIMF,  ATC  Sand-Clay,  ATC  DSTS).  The  complicated 
response  of  geologic  materials  to  dynamic  compression  is  of  fundamental  im¬ 
portance  to  understanding  the  energy  transmitted  to  the  surface  by  buried 
explosives.  We  measured  the  shock  states  induced  via  planar  impact  exper¬ 
iments  on  the  Harvard  40-mm  gas  gun.  Shock  wave  velocities  in  the  soil  sam¬ 
ples  were  measured  using  both  VISAR  and  piezoelectric  pins.  The  soils  were 
composed  primarily  of  quartz  with  different  mass  fractions  of  phyllosilicates 
and  amorphous  material.  Using  initial  particle  sizes  ranging  from  150  to  300 
microns,  the  samples  were  pressed  to  densities  ranging  from  1.7  to  1.93  g  cm“^ 
(about  25-30%  porous).  Water-saturated  samples  had  densities  ranging  from 
2.23  to  2.33  g  cm“^.  We  hud  that  the  dry  soils  have  a  linear  Ug  —  Up  rela¬ 
tion  that  is  similar  to  dry  quartz  sand  with  the  same  initial  density.  The  water- 
saturated  samples  are  less  compressible  and  have  much  greater  scatter  in  shock 
velocities.  The  VISAR  measurement  records  the  dispersion  around  the  mean 
shock  state  that  arises  from  reflections  between  grains,  and  we  will  compare 
the  VISAR  data  to  mesoscale  hydrocode  simulations  of  the  experiment.  These 
data  will  be  used  to  generate  more  accurate  rheological  models  for  hydrocode 
simulations  of  the  shock  response  of  heterogeneous  granular  materials  in  the 
low-pressure  regime  (<  10  GPa). 
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1.  Introduction 

Buried  explosives  pose  a  serious  threat  to  both  personnel  and  equipment.  The  output 
of  such  blasts  is  a  function  of  the  depth  of  burial,  initial  stress,  and  soil  conditions.  How¬ 
ever,  the  influence  of  the  physical  properties  (e.g.,  density  and  water  content)  and  dynamic 
response  of  soils  on  the  blast  output  is  not  well  understood  or  quantihed.  The  limited 
amount  of  data  available  for  soils  and  other  porous/granular  geologic  materials  demon¬ 
strate  that  the  dynamic  response  is  strongly  dependent  on  material  and  water  content. 
Knowledge  of  the  both  loading  and  unloading  behavior  is  necessary  to  determine  terminal 
velocities  of  soils  after  explosive  loading.  In  addition,  current  numerical  techniques  are 
unable  to  predict  a  priori  the  dynamic  response  of  soils  under  the  wide  range  of  conditions 
encountered  in  the  held.  This  work  falls  under  the  solid  mechanics  research  area  in  the 
U.S.  Army  Research  Office  by  providing  basic  shock  response  data  on  materials  identihed 
to  be  of  strategic  interest. 

2.  Experimental  Technique 

Using  standard  references  soils  provided  by  Army  Research  Laboratory  (VIMF,  ATC 
Sand-Clay,  ATC  DSTS),  we  conducted  shock  loading  experiments  on  compacted  dry  sam¬ 
ples  and  water-saturated  samples. 

2.1.  Material  properties 

Our  experiments  were  conducted  on  material  that  was  sifted  to  grain  sizes  ranging  from 
150  to  300  microns.  The  natural  soils  were  sifted  to  enable  the  shock  wave  to  traverse 
multiple  grains  during  the  dynamic  experiment.  The  grain  size  is  constrained  by  the 
duration  of  uniaxial  how  on  a  40-mm  gun. 

Table  1  presents  the  the  mineral  abundances  for  the  three  soils  measured  using  an  X-ray 
dihraction  technique.  Table  2  shows  the  particle  density  which  is  measured  using  helium 
pycnometry.  Information  about  the  particle  density  as  well  as  the  bulk  density  is  used  to 
calculate  the  porosity  of  the  soil  samples. 

2.2.  Sample  preparation 

Dry  samples  were  prepared  in  a  26  mm  diameter  capsule.  To  promote  good  mixing 
of  the  soil  constituents  and  hat  surfaces  the  samples  are  hrst  pressed  by  hand,  then 
lightly  hammered,  and  hnally  pressed  to  1  ton  up  to  three  times.  The  amplitude  of  the 
pre-experiment  press  mimics  precompression  of  soils  under  road  conditions.  The  multiple 
pressing  allowed  for  rotation  of  grains  and  a  more  homogeneous  pre-experiment  compacted 
state. 
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Wet  samples  were  prepared  in  a  28  mm  diameter  capsule  using  the  following  three-step 
process: 

1.  Dry  soil  was  added  to  the  capsule. 

2.  Water  was  added  to  the  capsule  in  excess  and  allowed  to  sit  for  a  period  of  approx¬ 
imately  20  hours. 

3.  A  specihcally  constructed  polycarbonate  piston  window  with  semi-cylindrical  cuts 
through  the  side  was  used  to  press  the  mixture  to  maximum  unshocked  density.  This 
setup  allowed  excess  water  to  escape  from  the  capsule  during  this  compression  process. 

The  water-saturated  samples  had  densities  ranging  from  2. 2-2. 6  g  cm“^,  corresponding 
to  the  25-30%  porosity  (Table  2). 

2.3.  Plate-impact  details 

Plate  impact  experiments  were  performed  on  the  40-mm  gas  gun  facility  at  Harvard 
University.  Experiments  measured  the  shock  properties  of  three  soils  (VIMF,  ATC  Sand- 
Clay,  and  ATC  DSTS)  under  both  dry  and  water-saturated  conditions.  For  each  soil 
sample  iteration,  we  conducted  impact  tests  at  approximately  300,  500,  800,  and  1200  m 
s-h 

Figure  1  shows  the  plate  conhguration  which  includes  the  soil  sample  mounted  between 
an  aluminum  2024  driver  and  a  PMMA  window.  The  soil  is  contained  within  a  poly¬ 
carbonate  capsule.  Shock  measurements  were  made  using  both  piezoelectric  pins  and  a 
velocity  interferometry  system  for  any  reflector  (VISAR)  [Barker  and  Hollenbach,  1972]. 
Both  the  piezoelectric  pins  and  the  VISAR  are  used  to  detect  shock  arrival  times  on  the 
aluminum  driver  and  the  soil  sample.  In  some  cases,  the  VISAR  is  used  to  measure  the 
particle  velocity  history  at  the  soil/PMMA  interface.  At  times,  the  VISAR  signal  de¬ 
graded  quickly,  in  some  cases  from  abrasion  of  the  vapor-deposited  aluminum  layer  on  the 
PMMA  window.  Thicker  aluminum  layers  were  used  more  successfully.  In  conjunction, 
these  two  techniques  allow  us  to  make  redundant  measurements  of  particle  velocity  which 
improves  experimental  uncertainties. 

2.4.  Further  Experimental  Considerations 

Our  experimental  analysis  relies  on  the  assumption  that  our  measurements  are  being 
made  under  uniaxial  compression.  However,  this  is  only  true  within  a  certain  spatial  and 
temporal  regime  in  the  experiment.  The  goal  of  this  analysis  will  be  to  show  that  with 
la  uncertainty,  all  measurements  are  made  under  uniaxial  compression. 

At  the  outer  diameter  of  the  soil  sample,  there  is  an  impedance  mismatch  between  the 
polycarbonate  capsule  and  the  soil  which  will  cause  lateral  release  waves  to  propogate 
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radially  inward  towards  the  center  of  the  sample.  The  leading  edge  of  these  waves  will 
propogate  at  the  bulk  sound  velocity  in  the  shock  state  of  the  soil,  given  by: 


cl=  I 


Cb  = 


(1) 


Using  a  Mie-Griineisen  equation  of  state,  we  can  use  the  known  Hugoniot  of  the  soils, 
calculated  from  the  impedance  match,  to  approximate  this  derivative  [Duffy  and  Ahrens, 
1992], 

'dP\ 

)  H 

Since  we  are  taking  our  measurements  at  the  back  surface  of  the  target,  we  use  expression 
for  the  diameter  on  the  rear  surface  where  the  compression  is  still  uniaxial,  derived  by 
[Swift  and  Kraus,  2008]: 

di  =  do  —  2/otan<h  (2) 


i-|i-U  ^2' 

Po  py  2  ^ 


+^V 

2p  / 


where  di  is  the  region  not  affected  by  lateral  release,  do  is  the  outer  diameter  of  the  target, 
and  lo  is  the  thickness  of  the  target,  and  <h  is  the  release  angle,  given  by: 


$  =  tan  ^ 


-  {Us-Upfy 


(3) 


Using  the  Hugoniot  of  the  soil  measured  by  the  impedance  match  calculation  and  a  Mie- 
Griineisen  parameter  of  1.2±0.3,  di  is  plotted  in  Figures  2-5  as  a  function  of  shock  velocity 
for  both  the  dry  and  wet  VIMF.  Error  bars  are  generated  using  the  Taylor  series  method 
discussed  in  Section  3.3.3. 


3.  Results 

3.1.  Shock  Hugoniot  States 

Given  an  initial  thermodynamic  state,  the  Rankine-Hugoniot  conservation  equations 
govern  the  locus  of  thermodynamic  states  that  are  attainable  by  a  single  shock  process. 
These  relations  require  the  conservation  of  mass,  momentum  and  energy  across  a  shock 
front  and  are  given  by: 


/  Ui\ 

Upi  UpO  —  Us  1  1  j 

(4) 

Pi  —  Pq  =  —{upi  —  Upo) 

ro 

(5) 

-Eo  =  ^{P,-Po){Vo-V,) 

(6) 

In  the  above  equations,  the  0  subscripts  denote  the  unshocked  state  and  the  1  subscripts 
denote  the  shocked  state.  In  our  experimental  setup,  the  initial  particle  velocity  is  zero 
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and  the  initial  pressnre  is  approximately  zero.  The  initial  specihc  volnme,  Vq  is  calculated 
by  careful  measurement  of  the  mass  and  volume  of  a  sample  element.  Thus,  to  fully 
constrain  the  problem  and  hnd  the  relationship  between  P,  V,  and  E  on  the  Hugoniot 
we  need  to  experimentally  measure  the  relationship  between  the  shock  velocity,  Us,  and 
the  particle  velocity,  Upi. 

3.1.1.  Impedance  Match  Method 

For  the  purposes  of  this  analysis,  parameters  relating  to  the  aluminum  driver,  the  soil 
sample,  and  the  PMMA  window  will  be  denoted  by  the  subscripts  a,  b,  and  c,  respectively. 
At  time  t=0,  the  aluminum  flyer  impacts  the  driver  generating  a  shock  wave  which  excites 
a  particle  velocity  of  —  ^  the  flyer  velocity  and  |  the  flyer  velocity  in  the  flyer  and  driver 
respectively.  When  the  shock  wave  reaches  the  driver/soil  interface,  there  is  a  partial 
release  of  pressure  along  the  driver’s  isentrope  due  to  the  lower  impedance  of  the  soil.  We 
assume  that  this  isentrope  is  well  approximated  by  the  principal  Hugoniot  of  aluminum 
reflected  about  the  axis  in  P  —  Up  space  where  the  particle  velocity  is  equal  to  |  the 
flyer  velocity.  This  is  generally  a  good  assumption  in  the  relatively  low  pressure  regime 
attained  in  these  experiments.  Using  both  piezoelectric  pins  and  a  velocity  interferometry 
system  for  any  reflector  (VISAR)  as  separate  non-correlated  measurements,  we  hnd  the 
time  at  which  the  shock  wave  reaches  the  driver/soil  interface.  Similarly,  we  measure 
the  time  at  which  the  shock  wave  reaches  the  soil/PMMA  interface.  This  allows  us  to 
measure  the  transit  time  of  the  shock  through  the  soil  sample.  By  measuring  the  sample 
thickness  prior  to  the  shock  experiment,  we  are  able  to  calculate  the  shock  velocity  in 
the  soil.  Shock  velocities  are  made  more  accurate  by  applying  a  flyer  plate  tilt  correction 
which  accounts  for  the  deviation  of  the  plane  of  the  driver  from  the  plane  of  the  projectile. 

The  corresponding  particle  velocity  is  obtained  using  two  separate  methods.  The  hrst 
is  an  impedance  matching  method,  which  requires  that  the  driver  release  to  the  state 
that  intersects  the  Rayleigh  line  of  the  soil  (Figure  6).  From  Equation  5,  the  Hugoniot  of 
aluminum  in  the  P  —  Up  plane  is: 

P  PaUsa'^pbl 

From  experimental  data  previously  collected  on  aluminum,  we  know  that  there  is  a  linear 
relationship  between  particle  velocity  and  shock  velocity  within  our  experimental  domain. 
Plugging  in  Ug  =  Ca  +  SaUp,  where  Ca  is  an  experimentally  measured  value  related  to  the 
bulk  sound  speed  of  aluminum,  and  Sa  is  an  experimentally  measured  constant,  we  can 
rewrite  our  expression  for  the  Hugoniot  of  aluminum; 

P  PaiCa  T  Sa'^pbl')l^pbl 
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If  we  assume  release  along  the  reflected  Hugoniot  of  aluminum,  then  the  expression  for 
the  release  path  is: 


P  —  Pa[Ca  "h  Sa(Vf  Ltpfel) 

The  impedance  match  constraint  tells  us  that: 

Pa[Ca  T  ^p)  PbU sb'^pbl 

Solving  for  Uphi-,  we  hnd  that: 


CaPa  +  ‘^SaVfPa  +  UsbPf,  +  a/  CaPa  +  ‘^^aUgbPaPb  +  ^iSaUgbV  f  PaPb  +  UsbPb 

^b 


(7) 


Aluminum  2024  was  chosen  as  a  flyer  and  driver  because  its  shock  properties  are  well 
characterized,  which  reduces  the  uncertainty  in  the  impedance  match  calculation  (by 
minimizing  the  uncertainty  in  the  driver  Hugoniot  parameters  Ca  and  Sa)-  Additionally, 
the  impedance  of  aluminum  is  a  better  match  to  the  impedance  of  our  samples  than  other 
common  impactors  such  as  steel,  which  minimizes  the  uncertainty  due  to  the  assumption 
that  the  isentrope  of  the  driver  is  well  approximated  by  the  Hugoniot.  We  use  the  data 
in  Marsh  [1980]  with  Up  <  900  m  s“^  to  £t  an  aluminum  2024  Hugoniot  in  the  pressure 
range  of  interest:  Us  =  5.362  +  1.3026Mp.  Using  the  impedance  match  method,  we  found 
the  following  Ug  —  Up  relations: 


^  0.5270  +  2.1545upfei 
0.3249  +  2.6445upfei 
0.3050  +  2.3226upfei 

2.0818  +  2.1464upfei 
2.1962  +  1.5135upfei 
,  2.3583  +  2.3608Mpfei 


Dry  VIMF 

Dry  ATC  Sand-Clay 

Dry  ATC  DSTS 

Wet  VIMF 

Wet  ATC  Sand-Clay 

Wet  ATC  DSTS 


where  the  £t  parameters  had  the  correlation  coefficients: 


(8) 


-0.960 

Dry 

VIMF 

-0.935 

Dry 

ATC  Sand- Clay 

-0.914 

Dry 

ATC  DSTS 

-0.938 

Wet 

VIMF 

-0.956 

Wet 

ATC  Sand-Clay 

-0.939 

Wet 

ATC  DSTS 

3.2.  VISAR  method 

Using  a  VISAR,  we  can  directly  measure  the  particle  velocity  at  the  soil  sample/  PMMA 
interface,  which  is  a  state  on  the  Hugoniot  of  PMMA  and  the  release  isentrope  of  the 
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soil.  To  find  the  particle  velocity  in  the  shock  state  of  the  soil,  we  need  to  know  both 
the  Hugoniot  of  PMMA  (shown  in  red  in  Figure  6)  as  well  as  an  expression  for  the 
isentropic  pressure  release  path  from  the  soil  to  the  PMMA  (shown  in  green  in  Figure  6). 
To  determine  the  Hugoniot  of  PMMA,  we  used  a  cubic  £t  to  the  dataset  published  by 
[Marsh,  1980]  at  relevant  particle  velocities  {up  <  1  km  s“^). 


Uc  =  Cc  +  ScUp  +  S'^up^  +  S”up^  (10) 

Where  Cc  =  2.8,  Sc  =  2.1,  S''  =  —2.9,  and  S'"  =  2.1.  In  P  —  Up  space,  the  general 
expression  for  the  Hugoniot  of  PMMA  can  be  written  as: 

Phc  =  PcUscUp  =  Pc{Cc  +  ScUp  +  S''mp^  +  S”up^)up  (11) 

Next,  we  assume  that  the  isentropic  release  is  well  approximated  by  the  reflected  Hugoniot 
of  the  soil.  We  use  the  Us  —  Up  relation  calculated  from  the  impedance  match  method  to 
make  this  correction.  This  means  that  the  particle  velocity  measured  from  the  VISAR 
method  is  not  completely  uncorrelated  from  the  impedance  match  method.  However, 
the  particle  velocity  correction  due  to  the  impedance  mismatch  between  the  soil  and  the 
PMMA  is  small  relative  to  the  total  particle  velocity,  so  we  treat  the  particle  velocity 
measured  from  the  VISAR  method  as  approximately  uncorrelated  to  the  particle  velocity 
measured  from  the  impedance  match  method.  From  the  impedance  match  data,  we  see 
that  the  Hugoniot  of  the  soil  is  well  fit  by  a  line,  so  the  Hugoniot  of  the  soil  is: 

Pnb  —  Pb{Ch  +  SfyUp^Up  (12) 

So  the  expression  for  the  approximate  isentrope  is: 

Psb  ~  Pb[Cb  T  >S'fe(2'Upft2  Ltp)](2'Upfe2  Lip)  (13) 

Where  Upb  is  the  axis  of  reflection,  which  is  equal  to  the  particle  velocity  in  the  shock 
state  of  the  soil.  By  the  impedance  match  criterion,  we  know  that  the  pressure  on  the 
isentrope  of  the  soil  must  be  equal  to  the  pressure  on  the  Hugoniot  of  the  PMMA: 

Phc  =  Psb  (14) 


Plugging  Equation  11  and  13  in  to  14  and  solving  for  Upb2i  we  get  the  following  expression 
for  the  particle  velocity  in  the  shock  state  of  the  soil  sample  as  a  function  of  the  measured 
particle  velocity  and  the  Hugoniot  parameters  of  the  soil  sample  and  PMMA: 


Upb2 


-CbPb  +  2SbPbUp  +  \/pb{Clpb  +  ASbUp{Cc  +  Up{Sc  +  Up{S'^  +  S'^Up)))pc) 

^SbPb 


Using  the  VISAR  technique,  we  calculated  the  following  Ug  —  Up  relations: 


(15) 
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us2  =  {  2.5863  +  1.6950Mpb2  Wet  ATC-DSTS  (16) 

with  correlation  coefficient: 

P2  =  {  -0.934  Wet  ATC-DSTS  (17) 

3.3.  Shock  Hugoniot  Uncertainties 

Due  to  the  limited  number  of  experiments  we  were  able  to  conduct  for  each  soil  sample, 
a  rigorous  quantihcation  of  experimental  uncertainties  was  done.  The  goals  of  this  analysis 
are  to: 

1.  Relate  the  uncertainties  in  the  measured  parameters  to  the  uncertainties  in  the 
calculated  Hugoniot  parameters  from  Equations  4,  5,  7,  and  15. 

2.  Demonstrate  that  the  VISAR  Hugoniot  parameters  are  approximately  independent 
of  the  impedance  match  Hugoniot  parameters. 

3.  Dehne  the  important  experimental  parameters,  to  serve  as  a  way  to  reduce  uncer¬ 
tainties  in  future  plate  impact  experiments. 

4.  Relate  the  uncertainties  in  the  experimental  Hugoniot  parameters  to  the  uncertain¬ 
ties  in  the  equation  of  state  variables  (the  coefficients  of  the  Ug  —  Up  relation). 

The  data  shown  in  Tables  3-8  are  presented  as  unique  measurements  with  1-dimensional 
uncorrelated  uncertainties.  In  reality,  the  uncertainties  in  each  calculated  parameter  are 
correlated  and  fall  an  an  ellipse  in  n-dimensional  space,  where  n  is  the  number  of  cal¬ 
culated  parameters.  The  following  analysis  presents  the  variance  and  covariance  of  our 
experimental  calculations. 

3.3.1.  Measurement  Uncertainties 

1.  Uncertainties  in  shock  arrival  times  are  taken  to  be  75%  of  the  signal  rise  time  on 
the  upper  end  and  25%  of  the  signal  rise  time  on  the  lower  end  for  both  the  VISAR  and 
the  piezoelectric  pins.  The  shock  breakout  time  is  taken  to  be  50%  of  the  signal  rise  time. 
Uncertainties  in  the  shock  arrival  times  are  generally  on  the  order  of  ~10  ns. 

2.  Uncertainties  in  the  sample  density  are  calculated  using  the  published  accuracies  of 
our  micrometer  and  scale.  Taking  several  repeated  measurement  helps  to  increase  the 
precision  of  our  density  calculation. 

3.  Uncertainties  in  the  flyer  velocity  are  taken  to  be  75%  and  25%  of  the  voltage  drop 
caused  by  the  flyer  interrupting  the  laser  beam  lines.  The  flyer  arrival  time  is  taken  to  be 
50%  of  the  voltage  drop. 

3.3.2.  Hugoniot  Uncertainties 
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Let  fn{xi,  he  a  set  of  k  linearly  independent  functions  of  i  variables  (n  =  1, /c). 

We  have  k  =  6  functions  corresponding  to  Equations  4,  5,  7,  15  and  i  =  14  unique 
variables.  Given  the  variance-covariance  matrix  of  the  independent  variables  denoted  Sx, 
the  goal  is  to  hnd  the  variance-covariance  matrix  of  the  dependent  variables  Sf ,  which  will 
give  us  both  the  uncertainties  in  our  calculated  parameters  as  well  as  the  strength  of  the 
correlation  between  the  values  calculated  from  the  impedance  match  and  those  calculated 
from  the  VISAR.  This  is  used  to  justify  presenting  each  method  as  a  unique  measurement. 
We  developed  two  approaches  to  the  problem.  First  we  will  use  an  analytical  solution 
using  a  hrst  order  Taylor  approximation,  which  allows  us  to  quantify  the  error  contribution 
from  each  measurement  to  the  total  uncertainty  in  the  Hugoniot  parameters.  Second  we 
employ  a  Monte  Carlo  technique  that  generates  a  numerical  solution  to  check  the  validity 
of  the  analytical  approach. 

3.3.3.  Taylor  Series  Method 
Consider  the  set  of  vector  valued  functions: 


fn  =  fn{Xi,X2,  ...,Xi)  (18) 

where  /„  are  the  set  of  experimental  results  determined  from  i  measured  variables  Xi.  Each 
of  the  measurements  has  a  total  uncertainty  which  is  equal  to  the  difference  between  the 
true  value  of  the  variable  and  the  measurement  value.  The  total  uncertainty  can  be 
decomposed  into  both  random  and  systematic  components,  denoted  in  the  following  way: 


measured 


Xi  +  SXi  + 


true  sys  error  rand  error 


(19) 


Next  we  approximate  the  set  of  functions  /„  by  Taylor  expansion: 


/n  —  /n  +  ~ 

We  can  rewrite  this  in  terms  of  the  random  and  systematic  errors: 

dfn 


dxi 


\x,{sxi  -f  exi) 


(20) 


There  are  two  fundamental  assumptions  associated  with  Equation  20.  The  first  assump¬ 
tion  is  that  the  higher  order  terms  are  negligible,  which  is  reasonable  when  the  measure¬ 
ment  error  is  small  since  the  next  term  in  the  series  expansion  contains  the  square  of  the 
measurement  errors.  This  assumption  is  is  verified  by  the  Monte  Carlo  analysis.  The 
second  assumption  is  that  the  derivative  of  /„  evaluated  at  7*  is  a  good  approximation  of 
the  derivative  of  /„  at  Xi,  which  cannot  be  rectihed  by  the  Monte  Carlo  analysis,  since  we 
dont  know  the  values  of  Xi. 
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We  are  interested  in  seeing  how  the  distribution  of  errors  in  Xi  maps  to  the  distribution 
of  errors  in  /„,  so  square  both  sides  of  the  equation  to  find  the  variance  of  the  distribution 
about  fn-  We  assume  that  the  systematic  errors  will  be  small  compared  to  the  random 
errors,  so  we  set  exi  to  0.  However,  from  this  analysis,  we  can  look  at  the  prefactor  of  each 
of  the  systematic  error  terms  to  understand  how  much  a  systematic  error  in  the  variable 
Xi  will  affect  the  uncertainty  in  /„. 

2,  (dU  V 


2  f  dfn  df\ 

In  =  -T^SXiSX 


dxi 


^  dxj 


(21) 


Where  is  the  transpose  of  This  gives  us  the  desired  result  in  index  notation. 


dxi 


Notice  that  is  just  the  Jacobian  of  /„,  so  we  can  write  this  result  in  matrix  form: 


(22) 

3.3.4.  Monte  Carlo  method 

Here  we  take  a  numerical  approach  to  the  same  problem.  In  addition  to  the  value  of  each 
measurement,  we  know  the  covariance  matrix  Ex,  as  discussed  in  Section  3.3.1.  We  use 
this  information  to  generate  random  points  from  a  multivariate  normal  distribution.  This 
gives  us  an  array  of  observations  of  /„,  denoted  fnk  (the  k*'^  observation  of  the  parameter 
fn-  From  this  set  of  points,  we  calculate  the  covariance  matrix  of  /„: 

1  ^ 

^  -  Ji)iU  -  h) 

i 

For  large  enough  n,  converges  to  the  true  value  of  Sf  provided  that  there  are  no 

systematic  errors  which  are  unaccounted  for  in  the  input  covariance  matrix  Ex. 

3.3.5.  Discussion 

We  have  conhrmed  that  the  covariance  matrix  calculated  from  the  Taylor  series  method 
is  in  good  agreement  with  the  covariance  matrix  calculated  from  the  Monte  Carlo  method. 
In  this  paper,  we  present  the  uncertainties  calculated  from  the  Taylor  series  approach, 
as  the  analytical  method  allowed  us  to  decompose  the  total  uncertainty  in  /„  into  its 
component  uncertainties  from  the  variables  Xi  which  helped  to  verify  that  the  propagated 
uncertainties  made  sense. 

We  see  that  the  parameters  calculated  from  the  impedance  match  technique  are  inde¬ 
pendent  of  those  from  the  VISAR  technique.  This  is  as  expected  since  the  impedance 
match  calculations  relied  only  on  independent  measurements.  However,  the  uncertainties 
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in  the  parameters  calculated  from  the  VISAR  technique  are  correlated  to  the  those  of 
the  impedance  match  technique  through  the  Hugoniot  parameters  Cg  and  Sg  calculated 
from  the  impedance  match  technique.  To  evaluate  the  strength  of  the  linear  correlation 
between  the  particle  velocity  from  the  impedance  match  {upi)  and  the  particle  velocity 
calculated  from  the  VISAR  technique  (MP2),  we  use  the  Pearson  product-moment  correla¬ 
tion  coefficient,  which  gives  values  between  —1  and  1,  where  a  value  of  0  is  no  correlation 
between  the  two  variables.  The  correlation  coefficient  is  dehned  as: 

cov(a;,  y) 


As  shown  in  the  example  calculation,  the  correlation  coefficient  between  Upbi  and  Upb2  is 
on  the  order  of  ~  0.001,  which  implies  that  the  variables  are  approximately  independent. 
This  justihes  treating  each  as  an  independent  measurement. 

3.3.6.  Hugoniot  Linear  Fit  Uncertainties 
The  linear  hts  to  the  Ug  —  Up  data  are  calculated  using  the  weighted  total  least  squares 
method  described  by  [Krystek  and  Anton,  2007],  which  hts  a  straight  line  to  data  with 
uncertainties  in  both  the  independent  {up)  and  dependent  {Ug)  variables.  This  method 
also  allows  us  to  generate  the  covariance  matrix  for  the  £t  parameters  shown  here. 


Dry  VIMF  1 
Dry  ATC-SandClay  1 
Dry  ATC-DSTS  1 
Wet  VIMF  1 
Wet  ATC-SandClay  1 
Wet  ATC-DSTS  1 
Wet  ATC-DSTS  2 


0.0006  -0.0085 
-0.0085  0.0014 

(23) 

0.0011  -0.0017' 

(24) 

-0.0017  0.0033 

0.0002  -0.0003' 

(25) 

-0.0003  0.0006 

0.0121  -0.0204' 

(26) 

-0.0204  0.0391 

0.0057  -0.0092' 

(27) 

-0.0092  0.0163 

0.0019  -0.0036' 

(28) 

-0.0036  0.0080 

0.0671  -0.1286' 

(29) 

-0.1286  0.0671 

(30) 

3.4.  Snmmary  of  Hugoniot  Measurements 

The  measured  Hugoniot  states  for  each  soil  are  presented  in  Figures  7-12  and  Tables  3-8. 
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4.  Discussion  and  future  work 

In  all  cases,  water-saturated  soils  have  a  larger  shock  velocity  and  higher  impedance 
compared  to  the  dry  samples.  The  wet  soil  Hugoniot  data  are  in  excellent  agreement 
with  the  ice-saturated  quartz  sand  data  from  Kraus  et  ah  [2010].  However,  we  found 
considerable  scatter  in  water-saturated  samples.  The  scatter  was  more  pronounced  in 
singly-pressed  samples  (primarily  ATC  Sand-Clay).  Subsequently,  multiply-pressed  sam¬ 
ples  signihcantly  reduced  the  scatter  bewteen  the  different  measurement  techniques  on 
a  single  sample.  In  one  experiment,  a  water-saturated  sample  had  the  water  pumped 
out  in  the  target  chamber  immediately  before  the  experiment.  This  experiment  had  a 
shock  velocity  nearly  identical  to  the  dry  sample.  A  survey  of  the  literature  hnds  consid¬ 
erable  scatter  in  wet  soil  or  sand  samples.  In  some  cases,  the  wet  shock  velocity  is  very 
close  to  the  dry  shock  velocity  (e.g..  Chapman  et  ah  [2006]),  and  we  suggest  that  water 
loss  may  have  occurred  in  these  cases.  We  added  a  capacitance  sensor  to  the  wet  target 
conhguration  to  check  for  water  loss  in  the  following  experiments. 

Detailed  simulations  of  the  wave  propagation  through  the  dry  and  wet  soils  for  compar¬ 
ison  to  the  VISAR  data  records  is  ongoing.  In  Figure  13,  the  particle  velocities  measured 
in  the  hnal  shock  state  display  a  characteristic  ’ringing’  pattern.  Based  on  our  ongoing 
simulations  using  the  CTH  shock  physics  code,  the  ringing  arises  from  wave  reflections 
between  the  grains.  Further  analyses  of  these  data  will  improve  the  modeling  of  hetero¬ 
geneous  granular  materials  (e.g.,  Vogler  et  al.  [2010]). 

5.  Presentations 

Kraus,  R.  G.,  M.  G.  Newman,  S.  T.  Stewart,  Hugoniot  Measurements  on  Heterogeneous 
Geologic  Materials,  Lunar  and  Planetary  Science  Gonference  43,  2680,  2012. 

Newman,  M.,  S.  T.  Stewart,  R.  G.  Kraus,  Hugoniot  Measurements  on  Dry  and  Water- 
Saturated  Soils,  American  Geophysical  Union  Fall  Meeting,  Abs.  PllA-1788,  2012. 

Journal  paper  in  preparation:  Newman,  M.  G.,  S.  T.  Stewart,  R.  G.  Kraus,  Hugoniot 
Measurements  on  Dry  and  Water-Saturated  Soils,  Journal  of  Geophysical  Research. 
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Table  1.  Mineral  abundances  (wt  %).  Quartz,  kaolinite,  muscovite  and  hematite  were 
identified  in  each  of  the  clays  using  x-ray  diffraction.  Amorphous  materials  may  include 
amorphous  silica,  Fe  oxyhydroxides,  and  poorly  crystalline  or  amorphous  clays. 


Soil 

Name 

Quartz 

Kaolinite 

Muscovite 

Hematite 

Amorphous 

1 

VIMF 

60.4 

13.2 

5.1 

0.3 

21.0 

2 

ATC  Sand-Clay 

67.2 

16.4 

7.2 

0.2 

9.0 

3 

ATC-DSTS 

56.5 

1.0 

11.0 

0.5 

31.0 

Table  2.  Soil  characteristics. 


Soil 

Name 

Crystal  density 
(g/cm^) 

Dry  bulk 
density  (g/cm^) 

Porosity 

Water-saturated 
bulk  density  (g/cm^) 

1 

VIMF 

2.539 

1.843 

0.274 

2.231 

2 

ATC  Sand-Clay 

2.566 

1.914 

0.254 

2.331 

3 

ATC-DSTS 

2.548 

1.750 

0.313 

2.297 
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Soil  Sample  (Dry  orWfet) 


Figure  1.  Target  schematic  (side  view)  for  determining  the  Hugoniot  state  via  measnrement 
of  the  shock  wave  transit  time.  The  VISAR  is  nsed  to  determine  if  any  corrections  need  to  be 
made  for  a  mnlti-wave  structure.  The  soil  experiments  utilized  simultaneous  piezoelectric  pins 
and  VISAR  at  2  to  5  points  on  the  sample. 


Figure  2.  The  diameter  over  which  the  experiment  remains  one  dimensional  with  1-a  error  as 
a  function  of  shock  velocity  for  the  dry  soil.  The  location  of  the  pins  and  the  visar  are  shown  in 
red  and  green  respectively. 
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Figure  3.  Release  angle  as  a  function  of  shock  velocity  with  l-cx  error  over  the  relevant 
experimental  domain  for  the  dry  soil. 


Figure  4.  The  diameter  over  which  the  experiment  remains  one  dimensional  with  1-cr  error  as 
a  function  of  shock  velocity  for  the  wet  soil.  The  location  of  the  pins  and  the  visar  are  shown  in 
red  and  green  respectively. 
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Figure  5.  Release  angle  as  a  function  of  shock  velocity  with  1-a  error  over  the  relevant 
experimental  domain  for  the  wet  soil. 


Figure  6.  Summary  of  Hugoniot  measurements  using  two  methods.  The  open  circle  represents 
the  value  calculated  by  the  impedance  match  method  and  the  closed  circle  represents  the  value 
measured  directly  by  the  VISAR. 
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Figure  7.  Principal  Hugoniot  of  the  dry  VIMF.  Points  shown  in  red  are  collected  from  the 
pin  impedance  match  and  points  shown  in  blue  are  collected  from  the  VISAR  impedance  match. 
Error  bars  are  l-cr  calculated  from  the  Taylor  series  method. 


Wet  VIMF 


Figure  8.  Principal  Hugoniot  of  the  wet  VIMF.  Points  shown  in  red  are  collected  from  the  pin 
impedance  match,  points  shown  in  blue  are  collected  from  the  VISAR  impedance  match,  and 
points  shown  in  black  are  collected  from  the  VISAR  wave  profile.  Error  bars  are  1-a  calculated 
from  the  Taylor  series  method.  The  points  shown  in  pink  and  light  blue  fall  on  the  dry  VIMF 
hugoniot,  indicating  that  the  sample  lost  water  before  the  shot. 
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Figure  9.  Principal  Hugoniot  of  the  dry  ATC  Sand-Clay. 


Wet  ATC-SandClay 


Figure  10.  Principal  Hngoniot  of  the  wet  ATC  Sand-Clay. 
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Figure  11.  Principal  Hugoniot  of  the  dry  ATC  DSTS.  The  systematic  lower  particle  velocity 
measured  from  the  VISAR  wave  prohle  (shown  as  black  squares)  compared  to  the  impedance 
match  (shown  in  blue  for  the  VISAR  and  red  for  the  pins),  suggests  stress  relaxation  after  the 
initial  shock. 


Wet  ATC-DSTS 


Figure  12.  Principal  Hugoniot  of  the  wet  ATC  DSTS.  The  £t  to  the  VISAR  data  is  shown  in 
black  and  the  £t  to  the  impedance  match  data  is  shown  in  blue. 
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Time  (microsec) 

Figure  13.  Examples  of  VISAR  wave  profiles  at  the  water-saturated  VIMF  soil-PMMA 
interface.  Impact  velocities  were  528  and  796  m/s.  The  ’ringing’  at  the  top  of  the  wave  arises 
from  pressure  variations  in  the  heterogeneous  sample. 


Table  3.  Dry  VIMF  Impedance  Match  Data 


Shot 

Measurement 

Us 

Up 

P 

V 

099 

VISAR 

2.51 

± 

.04 

0.93 

± 

.01 

4.23 

± 

.09 

0.348 

± 

.006 

J 

2.55 

± 

.02 

0.92 

± 

.01 

4.27 

± 

.07 

0.351 

± 

.004 

K 

2.50 

± 

.02 

0.93 

± 

.01 

4.21 

± 

.06 

0.347 

± 

.004 

L 

2.47 

± 

.03 

0.93 

± 

.01 

4.17 

± 

.07 

0.343 

± 

.004 

098 

VISAR 

2.07 

± 

.05 

0.710 

± 

.004 

2.68 

± 

.05 

0.360 

± 

.006 

I 

2.10 

± 

.02 

0.713 

± 

.007 

2.71 

± 

.04 

0.360 

± 

.003 

J 

2.10 

± 

.02 

0.713 

± 

.007 

2.71 

± 

.04 

0.360 

± 

.004 

K 

2.04 

± 

.02 

0.716 

± 

.007 

2.65 

± 

.04 

0.358 

± 

.004 

L 

2.06 

± 

.02 

0.715 

± 

.007 

2.67 

± 

.04 

0.360 

± 

.003 

101 

VISAR 

1.58 

± 

.05 

0.470 

± 

.003 

1.38 

± 

.04 

0.380 

± 

.006 

I 

1.54 

± 

.02 

0.471 

± 

.003 

1.35 

± 

.02 

0.373 

± 

.003 

J 

1.58 

± 

.02 

0.469 

± 

.003 

1.38 

± 

.02 

0.377 

± 

.003 

K 

1.54 

± 

.02 

0.471 

± 

.003 

1.35 

± 

.02 

0.373 

± 

.003 

L 

1.51 

± 

.02 

0.472 

± 

.003 

1.33 

± 

.02 

0.368 

± 

.003 

102 

VISAR 

1.03 

± 

.03 

0.270 

± 

.005 

0.52 

± 

.07 

0.39 

± 

.02 

I 

1.08 

± 

.04 

0.267 

± 

.002 

0.54 

± 

.02 

0.401 

± 

.006 

J 

1.10 

± 

.04 

0.266 

± 

.002 

0.55 

± 

.02 

0.404 

± 

.006 

L 

1.06 

± 

.06 

0.267 

± 

.002 

0.53 

± 

.03 

0.398 

± 

.009 
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Table  4.  Wet  VIMF  Impedance  Match  Data 


Shot 

Measurement 

Us 

Up 

P 

z/ 

110 

VISAR 

2.1C 

l± 

.2 

0.24 

± 

.01 

1.16 

± 

.11 

0.39 

± 

.03 

I 

2.82 

± 

.23 

0.226 

± 

.008 

1.44 

± 

.12 

0.41 

± 

.03 

J 

2.64 

± 

.18 

0.231 

± 

.007 

1.38 

± 

.11 

0.40 

± 

.03 

K 

2.44 

± 

.24 

0.236 

± 

.008 

1.30 

± 

.12 

0.40 

± 

.03 

L 

2.62 

± 

.16 

0.231 

± 

.007 

1.37 

± 

.10 

0.40 

± 

.03 

103 

VISAR 

3.68 

± 

.11 

0.79 

± 

.02 

6.43 

± 

.40 

0.35 

± 

.03 

I 

3.76 

± 

.09 

0.78 

± 

.03 

6.52 

± 

.42 

0.36 

± 

.03 

J 

3.84 

± 

.15 

0.78 

± 

.03 

6.62 

± 

.45 

0.36 

± 

.03 

K 

3.93 

± 

.15 

0.77 

± 

.03 

6.72 

± 

.45 

0.36 

± 

.03 

L 

3.95 

± 

.21 

0.77 

± 

.03 

6.75 

± 

.49 

0.36 

± 

.03 

105 

VISAR 

2.93 

± 

.18 

0.37 

± 

.01 

2.41 

± 

.19 

0.39 

± 

.03 

I 

3.10 

± 

.14 

0.37: 

±  . 

012 

2.51 

± 

.17 

0.40 

± 

.03 

J 

3.07 

± 

.15 

0.37: 

±  . 

012 

2.50 

± 

.18 

0.40 

± 

.03 

K 

2.95 

± 

.21 

0.37: 

±  . 

013 

2.43 

± 

.20 

0.40 

± 

.03 

L 

2.94 

± 

.17 

0.37: 

±  . 

012 

2.42 

± 

.18 

0.40 

± 

.03 

104 

VISAR 

2.63 

± 

.32 

0.58 

± 

.02 

3.37 

± 

.37 

0.35 

± 

.03 

J 

3.17 

± 

.10 

0.53 

± 

.02 

3.71 

± 

.25 

0.37 

± 

.03 

L 

3.00 

± 

.12 

0.54 

± 

.02 

3.57 

± 

.25 

0.37 

± 

.03 

106* 

VISAR 

0.92 

± 

.06 

0.28 

± 

.01 

0.56 

± 

.09 

0.31 

± 

.05 

I 

1.36 

± 

.05 

0.263 

± 

.009 

0.79 

± 

.12 

0.36 

± 

.06 

J 

1.38 

± 

.06 

0.262 

± 

.009 

0.80 

± 

.12 

0.37 

± 

.06 

K 

1.25 

± 

.08 

0.267 

± 

.009 

0.74 

± 

.11 

0.36 

± 

.06 

L 

1.26 

± 

.06 

0.266 

± 

.008 

0.75 

± 

.11 

0.36 

± 

.06 

Table  5.  Dry  ATC  Sand-Clay  Impedance  Match  Data 


Shot 

Measurement 

Us 

Up 

P 

V 

092 

VISAR 

1.11 

± 

.06 

0.234 

± 

.002 

0.50 

± 

.02 

0.409 

± 

.008 

Pin 

1.10 

± 

.04 

0.235 

± 

.002 

0.50 

± 

.02 

0.408 

± 

.006 

114 

VISAR 

1.34 

± 

.03 

0.463 

± 

.003 

1.07 

± 

.02 

0.379 

± 

.005 

I 

1.32 

± 

.02 

0.464 

± 

.003 

1.06 

± 

.02 

0.375 

± 

.004 

J 

1.33 

± 

.03 

0.463 

± 

.003 

1.07 

± 

.02 

0.377 

± 

.005 

K 

1.32 

± 

.03 

0.464 

± 

.003 

1.06 

± 

.02 

0.375 

± 

.005 

L 

1.39 

± 

.03 

0.460 

± 

.003 

1.11 

± 

.02 

0.388 

± 

.005 

091 

VISAR 

1.44 

± 

.02 

0.447 

± 

.003 

1.22 

± 

.02 

0.365 

± 

.003 

Pin 

1.45 

± 

.02 

0.446 

± 

.003 

1.22 

± 

.02 

0.367 

± 

.004 

088 

VISAR 

2.21 

± 

.04 

0.678 

± 

.008 

2.88 

± 

.06 

0.360 

± 

.005 

Pin 

2.27 

± 

.04 

0.674 

± 

.008 

2.95 

± 

.05 

0.365 

± 

.004 

087 

VISAR 

2.65 

± 

.05 

0.91 

± 

.01 

4.59 

± 

.09 

0.344 

± 

.005 

Pin 

2.73 

± 

.02 

0.90 

± 

.01 

4.70 

± 

.08 

0.350 

± 

.003 
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Table  6.  Wet  ATC  Sand-Clay  Impedance  Match  Data 


Shot 

Measurement 

Us 

Up 

P 

z/ 

109 

VISAR 

1.99 

± 

.13 

0.250 

± 

.005 

1.16 

± 

.06 

0.375 

± 

.005 

J 

3.03 

± 

.14 

0.223 

± 

.005 

1.57 

± 

.06 

0.397 

± 

.003 

K 

2.67 

± 

.08 

0.231 

± 

.004 

1.44 

± 

.04 

0.392 

± 

.002 

L 

2.604 

:± 

.07 

0.233 

± 

.004 

1.42 

± 

.03 

0.391 

± 

.002 

095 

VISAR 

2.30 

± 

.10 

0.410 

± 

.006 

2.08 

± 

.07 

0.373 

± 

.005 

J 

2.60 

± 

.13 

0.398 

± 

.007 

2.28 

± 

.09 

0.385 

± 

.005 

L 

3.17 

± 

.14 

0.376 

± 

.008 

2.62 

± 

.09 

0.400 

± 

.004 

097 

VISAR 

2.79 

± 

.15 

0.62 

± 

.01 

3.87 

± 

.16 

0.347 

± 

.007 

I 

3.16 

± 

.05 

0.60 

± 

.01 

4.22 

± 

.09 

0.361 

± 

.003 

J 

3.14 

± 

.10 

0.60 

± 

.01 

4.20 

± 

.12 

0.361 

± 

.004 

K 

3.18 

± 

.10 

0.59 

± 

.01 

4.24 

± 

.12 

0.362 

± 

.004 

L 

3.04 

± 

.09 

0.60 

± 

.01 

4.11 

± 

.11 

0.357 

± 

.004 

094 

VISAR 

3.12 

± 

.07 

0.84 

± 

.02 

5.86 

± 

.15 

0.326 

± 

.005 

J 

3.41 

± 

.08 

0.82 

± 

.02 

6.25 

± 

.17 

0.340 

± 

.004 

L 

3.64 

± 

.05 

0.80 

± 

.02 

6.53 

± 

.16 

0.348 

± 

.003 

Table  7.  Dry  ATC  DSTS  Impedance  Match  Data 


Shot 

Measurement 

Us 

Up 

P 

u 

119 

VISAR 

0.88 

± 

.11 

0.274 

± 

.003 

0.42 

± 

.05 

0.39 

± 

.02 

I 

0.89 

± 

.03 

0.273 

± 

.001 

0.43 

± 

.01 

0.40 

± 

.07 

J 

0.98 

± 

.01 

0.271 

± 

.001 

0.465 

± 

.006 

0.414 

± 

.003 

K 

0.95 

± 

.02 

0.273 

± 

.001 

0.43 

± 

.01 

0.397 

± 

.005 

L 

0.89 

± 

.03 

0.274 

± 

.001 

0.42 

± 

.01 

0.396 

± 

.006 

114 

VISAR 

1.32 

± 

.04 

0.434 

± 

.003 

1.06 

± 

.03 

0.375 

± 

.007 

I 

1.32 

± 

.02 

0.464 

± 

.003 

1.05 

± 

.01 

0.374 

± 

.003 

J 

1.33 

± 

.03 

0.463 

± 

.003 

1.07 

± 

.02 

0.377 

± 

.006 

K 

1.32 

± 

.03 

0.464 

± 

.003 

1.06 

± 

.02 

0.375 

± 

.005 

L 

1.39 

± 

.01 

0.461 

± 

.003 

1.11 

± 

.01 

0.386 

± 

.003 

121 

VISAR 

1.90 

± 

.03 

0.665 

± 

.006 

2.19 

± 

.04 

0.376 

± 

.005 

I 

1.80 

± 

.03 

0.671 

± 

.006 

2.09 

± 

.04 

0.363 

± 

.005 

J 

1.89 

± 

.02 

0.666 

± 

.006 

2.18 

± 

.03 

0.375 

± 

.003 

K 

1.87 

± 

.02 

0.667 

± 

.006 

2.15 

± 

.03 

0.371 

± 

.003 

L 

1.82 

± 

.04 

0.670 

± 

.006 

2.11 

± 

.05 

0.365 

± 

.006 

124 

VISAR 

2.47 

± 

.03 

0.94 

± 

.01 

4.16 

± 

.07 

0.346 

± 

.005 

I 

2.58 

± 

.01 

0.93 

± 

.01 

4.30 

± 

.06 

0.356 

± 

.003 

J 

2.44 

± 

.02 

0.94 

± 

.01 

4.12 

± 

.06 

0.343 

± 

.003 

K 

2.38 

± 

.02 

0.95 

± 

.01 

4.03 

± 

.06 

0.335 

± 

.004 

L 

2.53 

± 

.01 

0.93 

± 

.01 

4.24 

± 

.06 

0.352 

± 

.003 
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Table  8.  Wet  ATC  DSTS  Impedance  Match  Data 


Shot 

Measurement 

Us 

Up 

P 

u 

128 

VISAR 

2.97 

± 

.21 

0.223 

± 

.005 

1.51 

± 

.08 

0.406 

± 

.004 

I 

2.92 

± 

.06 

0.224 

± 

.003 

1.49 

± 

.03 

0.405 

± 

.002 

J 

2.92 

± 

.08 

0.224 

± 

.003 

1.49 

± 

.03 

0.405 

± 

.002 

K 

2.94 

± 

.08 

0.224 

± 

.003 

1.50 

± 

.03 

0.406 

± 

.002 

L 

2.97 

± 

.06 

0.223 

± 

.003 

1.51 

± 

.03 

0.406 

± 

.002 

127 

VISAR 

3.16 

± 

.21 

0.343 

± 

.009 

2.44 

± 

.12 

0.396 

± 

.005 

I 

3.12 

± 

.08 

0.345 

± 

.006 

2.42 

± 

.06 

0.395 

± 

.003 

J 

3.18 

± 

.07 

0.343 

± 

.006 

2.45 

± 

.06 

0.396 

± 

.002 

K 

3.14 

± 

.05 

0.344 

± 

.005 

2.43 

± 

.05 

0.396 

± 

.002 

L 

3.03 

± 

.04 

0.347 

± 

.005 

2.37 

± 

.04 

0.393 

± 

.002 

130 

VISAR 

3.43 

± 

.13 

0.536 

± 

.012 

4.20 

± 

.14 

0.367 

± 

.004 

I 

3.62 

± 

.04 

0.53 

± 

.01 

4.37 

± 

.09 

0.372 

± 

.002 

J 

3.68 

± 

.04 

0.52 

± 

.01 

4.42 

± 

.10 

0.373 

± 

.002 

K 

3.66 

± 

.04 

0.52 

± 

.01 

4.40 

± 

.09 

0.373 

± 

.002 

L 

3.58 

± 

.05 

0.53 

± 

.01 

4.34 

± 

.09 

0.371 

± 

.001 

129 

VISAR 

3.90 

± 

.16 

0.77 

± 

.02 

6.93 

± 

.25 

0.348 

± 

.005 

I 

4.14 

± 

.07 

0.75 

± 

.02 

7.21 

± 

.20 

0.354 

± 

.003 

J 

4.21 

± 

.08 

0.75 

± 

.02 

7.29 

± 

.21 

0.356 

± 

.003 

K 

4.19 

± 

.08 

0.75 

± 

.02 

7.27 

± 

.20 

0.355 

± 

.003 

L 

4.06 

± 

.05 

0.76 

± 

.02 

7.12 

± 

.18 

0.352 

± 

.002 
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